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INTRODUCTION 


Internal  blast  is  the  concussion  effect  in  a confined 
explosion.  Confined  air  is  suddenly  warmed,  generating  a 
pressure  rise  that  can  be  quite  damaging.  Examples  of  inter- 
nal blast  include  blast  from  explosions  in  ship  compartments, 
dust  explosions  in  grain  elevators  or  in  coal  mines,  explo- 
sions of  gasoline  fumes  in  motorboat  hulls  or  in  refinery 
tanks,  and  natural  gas  explosions  in  kitchens  or  in  entire 
buildings.  Internal  blast  may  also  account  for  a situation 
such  as  the  mysterious  disappearance  of  an  oil  tanker  at  sea. 


Suddenness  of  the  pressure  rise  in  these  internal  explo- 
sions results  from  a well-known  effect  of  pressure,  which  in- 
creases the  rate  of  chemical  reaction — that  of  combustion.  A 
small  pressure  rise  from  an  initial  temperature  rise  increases 
the  rate  of  combustion,  which  in  turn  increases  the  rate  of 
pressure  rise,  and  so  on,  until  some  limiting  condition  is 
reached.  For  internal  explosions,  the  limiting  condition  is 
the  peak  overpressure  for  internal  blast. 

With  internal  blast,  and  considering  both  conventional 
and  explosive  fuels  composed  of  carbon,  hydrogen,  nitrogen, 
and  oxygen,  and  including  argon  from  the  air,  there  are  five 
types  of  atoms  in  the  products.  These  products  may  involve 
as  many  as  12  different  chemical  species  including  unstable 
ones  such  as  hydroxyl  and  monatomic  hydrogen,  monatomic  oxy- 
gen, and  monatomic  nitrogen.  To  characterize  these  12  spe- 
cies, 12  independent  relations  are  required.  Five  of  these 
are  provided  by  the  law  of  conservation  of  mass,  one  for  each 
type  of  atom.  The  seven  additional  relations  are  obtained 
through  the  second  law  of  thermodynamics  where  the  minimum 
Gibbs  free  energy  restriction  is  utilized  in  the  form  of  chem- 
ical equilibrium  constants.  These  form  12  simultaneous  non- 
linear algebraic  equations  that  can  be  solved  by  a simple 
heuristic  method  that  is  described  and  illustrated,  and  for 
which  a small  tabletop  computer  is  quite  suitable. 


To  determine  blast  overpressure  in  an  internal  explosion, 
it  is  necessary  to  know  the  explosion  temperature  and  compo- 
sition of  products  at  that  temperature.  To  determine  explo- 
sion temperature,  the  conservation  of  energy  principle  is 
utilized.  The  energies  involved  must  include  both  thermal 
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energies  associated  with  chemical  composition.  When  these 
are  known,  the  blast  overpressure  is  computed  through  the 
ideal  gas  law.  Such  computations  are  all  described  by  ex- 
ample. 

Overpressures  computed  for  internal  explosions  with  TNT 
agree  well  with  experimental  values.  Such  theoretical  calcu- 
lations also  yield  items  not  readily  available  experimentally; 
e.g.,  explosion  flame  temperature  and  the  composition,  for- 
mula mass,  and  heat  capacity  of  explosion  products.  The  pro- 
gram also  allows  for  the  varying  effects  of  (1)  injecting  a 
fuel  at  high  velocity,  (2)  altitude  on  blast  overpressure, 
and  (3)  heat  transfer  from  a hot  flame  to  cool  confining 
walls . 


Overpressures  in  internal  explosions  depend  on  both  the 
amount  and  type  of  fuel  causing  the  explosion.  A fuel  can  be 
characterized  by  its  internal  blast  yield,  which  is  the  in- 
verse ratio  of  the  amount  of  a fuel  of  concern  to  the  amount 
of  a reference  fuel  that  provides  the  same  overpressure.  In- 
ternal blast  yields  have  been  computed  and  are  reported  for 
conventional  fuels  such  as  jet  fuel,  gases  such  as  ethylene 
and  methane,  and  explosive  fuels  such  as  bexanitrobenzene  and 
pentaerythritoltetranitrate. 


CHARACTERISTICS  OF  INTERNAL  BLAST 
PEAK  OVERPRESSURES 

Blast  in  an  internal  explosion  is  characterized  largely 
by  its  peak  overpressure  and  duration.  Peak  overpressures 
for  an  internal  explosion  ordinarily  are  only  a few  bars, 
perhaps  10  or  less,  and  are  low  compared  with  overpressures 
for  explosions  in  the  unconfined  atmosphere.  Maximum  peak 
overpressures  in  internal  explosions  occur  for  the  special 
case  of  (1)  no  venting  effects,  and  for  (2)  adiabatic  condi- 
tions with  no  cooling  effects  by  heat  transfer  to  confining 
walls.  Such  maximum  peak  overpressures  are  a primary  con- 
cern in  this  report . 

After  the  peak  overpressure  is  attained,  a quasi-expo- 
nential pressure-decay  phase  is  caused  by  cooling  effects  of 
confining  walls  and  by  gas  leakage  by  venting.  This  decay 
effect  is  relatively  slow  compared  with  rates  that  ordinarily 
occur  in  explosions  in  the  unconfined  atmosphere.  The  decay, 
however,  causes  the  pressure  to  approach  an  equilibrium  value 
that  for  vented  explosions  is  also  the  ambJent.  For  non- 
vented  explosions,  this  equilibrium  pressure  can  be  greater 
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than  ambient  if  additional  gas  volume  is  given  by  formation 
of  species  such  as  carbon  monoxide,  or  can  be  less  than  ambi- 
ent pressure  with  gas  diminishing  reactions  such  as  formation 
of  aluminum  oxide. 

The  slow  pressure-decay  rates  for  internal  explosions 
lead  to  relatively  long  overpressure  duration  times  that  are 
perhaps  as  long  as  a major  fraction  of  a second.  The  time 
integral  of  the  blast  overpressure  is  the  blast  impulse  per 
unit  area.  Blast  impulse  can  be  an  important  factor  in  blast 
damage  capability,  but  for  internal  blast  all  this  impulse 
may  not  be  effective  for  structural  damage.  This  situation 
arises  because  such  decay  time  may  exceed  a critical  time  for 
the  affected  structure.1 


EXPLOSIVE  VERSUS  CONVENTIONAL  FUELS 


Internal  blast  from  explosions  with  conventional  fuels 
is  the  direct  result  of  a simple  combustion  process.  Explo- 
sive fuels,  however,  differ  somewhat  in  that  an  initial  det- 
onation sets  up  a transient  explosive  shock,  one  that  is 
quickly  dissipated  and  can  often  be  neglected  in  the  overall 
energy  effects.  For  explosive  fuels  that  also  are  oxygen- 
deficient  (trinitrotoluene  (TNT),  for  example),  the  initial 
detonation  is  followed  by  combustion  in  an  afterburn;  it  is 
the  combined  detonation-combustion  reaction  that  generates 
internal  blast.  For  explosive  fuels  that  are  oxygen-rich 
(nitroglycerin,  for  example),  there  is  no  afterburn;  internal 
blast  is  only  result  of  the  warming  effect  of  a detonation. 

With  con  er.tional  fuels  the  internal  blast  overpressure 
reaches  a maximum  at  a fuel-air  ratio  that  corresponds  to 
optimum  utilization  of  the  oxygen  present.  Less  fuel  leaves 
unused  oxygen,  while  more  fuel  does  not  give  any  increased 
energy  release.  Excess  fuel  serves  only  to  increase  the 
amount  of  products  to  be  warmed,  and  the  reduced  temperature 
rise  that  results  gives  a decrease  in  developed  overpressure. 
The  characteristic  maximum  overpressure  for  conventional  fuels 
contrasts  with  the  behavior  of  explosive  fuels  that  carry 
their  own  oxygen;  developed  overpressures  for  explosive  fuels 
increase  monotonically  with  the  amount  of  fuel. 


1 R.  G.  S.  Sewell  and  G.  F.  Kinney.  "Response  of  Struc- 
tures to  Blast:  A New  Criterion,"  Ann.  .7.  jf.  Aoad.  Sai.t 
Vo  1 . 152  (1968),  p.  532. 
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THEORETICAL  BACKGROUND 


Efforts  to  describe  internal  blast  mathematically  date 
back  at  least  to  1945,  when  the  Los  Alamos  Scientific  Labo- 
ratory suggested  an  empirical  equation,2 

tP  - 3000(V//)  (1) 

in  then  acceptable  units  of  overpressure  hP  in  pounds  per 
square  inch,  mass  of  explosive  V in  pounds,  and  confined 
volume  of  air  V in  cubic  feet.  In  modern  metric  units  of 
bars,  kilograms,  and  cubic  metres  (steres)  this  can  be 
written  as 


up  = I3(v/y)  (2) 

The  metric  pressure  unit  is  the  bar,  defined  in  terms  of  SI 
units  as  10=  pascals  (newtons  per  square  metre),  or  as  one- 
tenth  of  a megapascal.  One  bar  corresponds  to  a pressure  of 
14.5  psi,  approximately  that  of  the  ordinary  atmosphere.  The 
units  of  the  term  (W/V)  of  the  equation,  kilograms  per  cubic 
metre,  are  almost  identical  with  ounces  per  cubic  foot. 

Overpressures  computed  through  empirical  Equation  2 for 
internal  explosion  with  TNT  in  air  are  shown  in  Figure  1. 

Also  shown  are  experimental  values  for  ordinary  sea  level 
conditions.3  Rough  agreement  is  indicated. 

Subsequently  ?n  equation  for  internal  blast  overpressure 
was  derived  using  the  ideal  gas  law,  but  ignoring  many  com- 
plicating effects.4  This  equation  states  that  for  adiabatic 
nonvented  conditions  the  overpressure,  pounds  per  square  inch, 
is 

CP  = 8.8 (H/V)  (3) 


’ Private  communication  between  James  Weeks  and  R.  W. 
Carlson,  Confinement  of  an  Explosion  by  a Steel  Vessel,  Los 
Alamos  Scientific  Report  LA-390  dated  1945. 

3 Hans  R'.  W.  Weibull.  "Pressures  Recorded  in  Partially 
Closed  Chambers  at  Explosions  of  TNT  Charges,"  Ann.  S.  'I. 
Aoad.  Sai.,  Vol.  152  (1968),  p.  357. 

4 National  Defense  Research  Committee.  Effects  of  Im- 
pact and  Explosion.  Vol.  1.  Washington,  D.  C.,  NDRC , 1946. 
P.  91.  (Summary  Report  of  Technical  Division  2,  Vol.  1, 
publication  UNCLASSIFIED.) 
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0 0.$  1.0  1.5 

TNT -AIR  RATIO,  kfl/m3 

FIGURE  1.  Internal  Blast  Overpressures 
lor  TNT-Air  Explosions. 

where  the  term  H/v  is  heat  of  combustion  in  kilocalories  per 

cubic  foot  of  product  gases.  In  metric  units  of  the  bar  and 

the  joule  per  cubic  metre,  this  equation  becomes 

AP  - 10-5 ( k - 1 )(£//)  (4) 

where  k is  the  ratio  of  heat  capacities  for  the  confined  gas, 

taken  as  1.4,  as  for  air  at  room  temperature.  Derivation  of 
this  equation  ignores  the  effects  of  gases  formed  in  the  ex- 
plosion, of  temperature  on  heat  capacity,  and  more  import- 
antly, the  energy  absorbing  effects  of  chemical  dissociations 
in  the  explosion-combustion  flame.  Equation  4 may  thus 
greatly  overstate  internal  blast  overpressures,  particularly 
for  explosions  with  appreciable  fuel-air  ratios  (Figure  1). 

An  alternative  approach  to  internal  blast  is  afforded 
by  thermodynamics.  Results  of  simplified  analysis  conform 
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approximately  to  experimental  values,  presumably  by  compensa- 
tion effects.5  An  improved  and  rigorous  thermodynamic  anal- 
ysis of  internal  blast  takes  into  account  the  complicating 
effects  described  above.  Particularly  important  among  these 
are  chemical  dissociations  that  set  an  effective  ceiling  on 
both  explosion  temperature  and  blast  overpressure . Combustion 
products  at  mechanical,  thermal  and  chemical  equilibrium  are 
considered  in  the  analysis,  omitting  transient  effects  of 
initial  explosive  shock. 

These  thermodynamic  calculations  are  based  on  an  ideal 
solution  of  ideal  gase^.  The  principles  of  conservation  of 
mass,  conservation  of  --nergy,  and  of  minimum  Gibbs  free  energy 
are  used.  These  considerations  are  taken  up  separately.  Rep- 
resentative results  from  such  calculations  are  shown  in  figure 
2.  In  this  log-log  plot,  computed  versus  experimental  over- 
pressures generated  for  TNT-air  internal  explosions  are  com- 
pared. Good  agreement  is  obtained  over  about  three  decades 
in  fuel-air  ratios,  a result  that  leads  to  confidence  in  the 
computations . 

In  thermodynamic  analysis  for  internal  blast,  a products 
composition  is  first  found  that  corresponds  to  some  more  or 
less  arbitrarily  selected  temperature.  Products  energy,  which 
must  include  both  thermal  and  chemical  aspects,  is  then  found 
for  each  component;  the  weighted  sum  of  these  gives  a total 
energy  for  the  products  system.  For  adiabatic  combustions 
this  total  energy  equals  that  for  the  original  fuel  and  air. 
Successive  trial  temperatures  are  then  selected  until  the 
adiabatic  condition  of  zero  energy  difference  is  found. 

Through  this  method  energy  accounting  proceeds  indirectly 
through  formation  energies  for  individual  components  rather 
than  through  direct  considerations  of  the  chemical  reactions 
themselves . 

Conventional  fuels  such  as  fuel  oil  contain  primarily 
carbon  and  hydrogen,  and  explosive  fuels  such  as  TNT  gener- 
ally contain  nitrogen  and  oxygen.  These  are  all  identified 
as  C-H-N-0  fuels.  Air  for  the  combustion  reaction  contains 


5 Air  Force  Systems  Command  for  Joint  Technical  Coor- 
dinating Group  for  Munitions  Effectiveness.  Internal  Blast 
Damage  Mechanisms  Computer  Program , by  James  F.  Proctor, 
Naval  Ordnance  Laboratory,  Silver  Spring,  Md.  Dayton,  Ohio, 
AFSC , ffright-Patterson  Air  Force  Base,  10  April  1973.  (Report 
61-JTCG/ME-73-3 , publication  UNCLASSIFIED.) 
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FIGURE  2.  Computed  (solid  line)  Versus  Experimental 
(points)  Overpressures  for  TNT-Air  Internal  Explosions. 


nitrogen,  oxygen,  and  argon,  so  that  five  different  atomic 
species  can  be  in  the  resulting  products  mixture. 

It  has  been  found  that  for  explosion  of  C-H-N-0  fuels 
in  air  there  may  be  as  many  as  12  different  chemical  species 
in  the  hot  gaseous  products . These  include  the  usual  com- 
bustion products:  carbon  dioxide,  CO2;  carbon  monoxide,  CO; 
water  vapor,  H2O ; hydrogen,  H2;  nitrogen,  N2;  oxygen,  O2; 
and  argon,  Ar.  Nitric  oxide  NO,  the  unstable  monatomics 
hydrogen  H,  oxygen  0,  and  nitrogen  N,  and  the  unstable  com- 
pound hydroxyl  OH  must  also  be  included.  These  12  species 
in  the  products  constitute  12  unknowns  so  that  12  independ- 
ent relations  are  required  to  find  the  amount  of  each  com- 
ponent present.  These  12  relations  are  given  in  the  fol- 
lowing discussions. 
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CONSERVATION  OF  MASS  RELATIONS 


The  conservation  of  mass  principle  Is  utilized  in  terms 
of  individual  elements.  Thus  all  carbon  in  the  products 
system  is  supplied  by  the  fuel  and  can  be  expressed  as  the 
number  of  moles  of  carbon  atoms,  eg.  This  equals  the  sum  of 
the  number  of  moles  of  carbon  dioxide,  CO2,  and  carbon  mon- 
oxide, CO,  in  the  products;  i.e..  » C02  + CO,  where  a 

chemical  formula  is  used  to  indicate  the  number  of  moles  of 
a designated  species.  This  type  of  equation  is  identified 
as  a material  balance.  Similarly  the  material  balance  for 
hydrogen  atoms,  ag,  equals  twice  the  number  of  moles  of 
water  vapor,  plus  twice  that  for  molecular  hydrogen,  plus 
that  for  monatomic  hydrogen  and  hydroxyl  OH.  Such  material 
balances  for  the  five  chemical  elements  are  written  as  fol- 
lows : 


Carbon 

aC  “ 

co2  + co 

(5) 

Hydrogen 

aE  - 

2H20  + 2H2  + H + OH 

(6) 

Oxygen 

a0  " 

2C02  + CO  + H20  + 20p 

+ 0 + OH  + NO  (7) 

Nitrogen 

aN  - 

2N2  + N + NO 

(8) 

Argon 

aAR  “ 

Ar 

(9) 

These  material  balances  represent  five  of  the  required 
12  relations;  the  remaining  seven  are  obtained  from  free 
energy  consideration. 


NOMINAL  PRODUCTS 


Preliminary  to  computation  of  the  composition  of  the 
products  in  an  internal  explosion  is  the  assignment  of  nom- 
inal values  for  the  number  of  moles  of  each  to  its  components 
There  are  several  different  sets  of  such  nominal  values,  all 
based  on  material  balance  relations.  In  the  method  selected, 
oxygen  is  assigned  first  to  conversion  of  carbon  into  carbon 
monoxide,  then  to  conversion  of  carbon  monoxide  into  carbon 
dioxide,  then  to  conversion  of  hydrogen  into  water,  and  if 
any  remains,  to  molecular  oxygen.  Such  nominal  compositions 
correspond  well  to  actual  combustions  at  lower  temperatures; 
at  high  temperatures  a more  realistic  nominal  composition  is 
obtained  by  considering  the  formation  of  water  before  the 
formation  of  carbon  dioxide. 
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STOICHIOMETRIC  REFERENCE  POINTS 

Various  stoichiometric  reference  points  can  be  estab- 
lished  in  terms  of  the  nominal  products  compositions.  The 
CO2-H2O  point  corresponds  to  the  fuol-air  ratio  for  complete 
combustion  to  CO2  and  H2O  without  excess  oxygen.  The  CO2- 
H2  point  represents  a lesser  oxygen  supply  (or  a greater 
fuel  supply)  where  the  nominal  products  contain  elemental 
hydrogen  rather  than  water  vapor.  The  CO-HoO  and  the  CO-H2 
points  are  similarly  defined.  Fuel-air  ratios  for  these 
four  stoichiometric  points  for  internal  explosions  with  TNT 
are  indicated  in  Figure  2.  Since  there  are  no  discontinui- 
ties in  the  pressure  curve  at  these  points,  they  do  not 
represent  situations  for  a sudden  change  in  products  compo- 
sition. A representative  computation  of  a stoichiometric 
reference  point  is  presented  in  Appendix  A. 

A fuel  supply  less  than  that  for  the  conventional  C02- 
H20  stoichiometric  point  corresponds  to  an  oxygen-rich  or 
fuel-lean  combustion;  a greater  fuel  supply  corresponds  to 
fuel-rich  conditions.  A fuel  supply  even  greater  than  that 
for  the  CO-H2  stoichiometric  point  may  also  correspond  to 
the  presence  of  solid  carbon  in  nominal  products.  The  CO-H2 
point  is  also  the  nominal  smoke  point;  a greater  fuel  supply 
can  be  expected  to  generate  a smoky  flame. 


GIBBS  FREE  ENERGY  RELATIONS 


The  second  law  of  thermodynamics  states  that  a closed 
system  in  equilibrium  at  uniform  pressure  and  temperature 
shows  a minimum  in  its  total  Gibbs  free  energy.  The  same 
rules  apply  to  the  quasistatic  situations  with  internal  ex- 
plosions. Gibbs  free  energy  values  for  particular  species 
are  conveniently  expressed  in  alternative  form  as  equilib- 
rium constants  of  formation;  data  are  available  in  sources 
such  as  the  JANAF  therraochemical  tables.6  These  tabulated 
data  are  readily  converted  into  conventional  thermodynamic 
equilibrium  corstants,  and  in  turn  into  working  equilibrium 
constants  in  terms  of  mole  numbers,  Kn  (Appendix  B). 


6 National  Bureau  of  Standards,  J All  AS  Thavnoohemical 
Tables , D,  R.  Stull  and  H.  Prophet,  Project  Directors. 
Washington,  D.  C. , NBS , 1971. 
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The  seven  mole-number  relations  required  here  and  ob- 
tained by  Gibbs  free  energy  considerations  are  shown  bslow. 
In  the  chemical  equation  to  the  left,  a chemical  formula 
identifies  a specific  component  in  a mixture  at  chemical 
equilibrium.  In  the  equilibrium  constant  equations  to  the 
right,  a chemical  formula  represents  the  number  of  moles  of 
that  component  present  in  the  equilibrium  mixture. 


Dissociation  of  carbon  dioxide: 

C02  - CO  + *j02  Kn  » CO  x 02i/C02  (10) 

Dissociation  of  water  vapor: 

H20  ■ H2  + ij02  Kn  - H2  x 02i/C02  (11) 

Formation  of  monatomic  hydrogen: 

ijH2  - H Hn  - H/H2  i (12) 

Formation  of  monatomic  oxygen: 

*l02  - 0 Kn  - 0/02i  (13) 

Formation  of  hydroxyl: 

W2  + h02  * OH  Kn  - OH/(d24  x H2*)  (14) 

Formation  of  nitric  oxide:  . 

*jN2  + h02  « NO  Kn  « N0/(nj*  X 02*J  (15) 

Formation  of  monatomic  nitrogen: 

4N2  - N Kn  - N/N2*  (16) 


Numerical  values  for  a mole  number  equilibrium  constant 
Xn  depend  on  the  species  involved,  the  temperature,  and  the 
volume  cf  the  system.  In  principle  such  extensive  equilib- 
rium data  could  be  stored  in  the  memory  of  a large  computer 
and  utilized  as  required.  The  needed  tabulations  can  be 
considered  condensed  into  relatively  simple  algebraic  equa- 
tions for  the  (logarithm  of  the)  equilibrium  constant  of 
formation  as  a function  of  absolute  temperature  T.  Thus 

log  Kf  - A + 3/Cr+C)  (17) 

where  constants  A and  3 are  essentially  those  of  the  classic 
van 1 1 Hoff  equation,  and  constant  C,  a curvature  constant, 
provides  for  small  discrepancies  l Jtween  experimental  and 
van't  Hoff  values.  Constants  of  Equation  17  are  tabulated 
in  Appendix  C for  the  seven  chemical  equilibria  in  the  in- 
ternal blast  generated  by  C-H-N-0  fuels.  The  temperature 
range  for  these  is  from  4000  K,  well  above  the  temperatures 
of  ordinary  explosions,  down  to  1500  K,  a quench  temperature 
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below  which  rearrangement  reactions  for  explosion  products 
are  so  sluggish  that  chemical  equilibrium  may  not  be  attained. 
The  general  suitability  of  these  equations  is  readily  veri- 
fied. 


DETERMINATION  OF  EQUILIBRIUM  COMPOSITIONS 


Twelve  simultaneous  equations  for  defining  the  equilib- 
rium compositions  described  above  (five  linear  material  bal- 
ances and  seven  nonlinear  equilibrium  expressions)  are  to  be 
solved.  Many  methods  for  solution,  some  of  them  quite  ele- 
gant, have  been  suggested.  A relatively  simple  heuristic 
method  of  successive  approximation,  where  each  approximation 
seeks  to  improve  on  previous  ones,  is  utilized  here.  The 
advantage  of  this  method  is  that  it  can  be  adapted  to  manual 
solution  and  is  also  readily  programmed  for  a small  computer. 

A suitable  algorithm  is  shown  in  Appendix  D. 

The  iteration  calculation  begins  with  the  nominal  com- 
position obtained  through  material  balance  considerations  as 
described  above.  Then  the  program  computes,  through  equilib- 
rium constant  expressions,  the  amount  of  other  species  (the 
minor  components)  that  would  be  in  chemical  equilibrium  with 
these  nominal  species  (the  major  components)  at  the  assigned 
temperature.  The  amounts  of  major  components  are  corrected 
for  minor  components  by  using  the  material  balance  equations. 
Ttte  minor  components  are  then  recomputed,  the  major  compon- 
ents recorrected,  and  so  on,  until  all  requirements  for  equi- 
librium have  been  met  as  shown  by  identical  successive  results 
in  the  computation. 

As  indicated  in  the  algorithm  of  Appendix  D,  details  for 
the  iteration  depend  on  whether  the  products  mixture  is  oxygen- 
rich  (the  nominal  products  contain  molecular  oxygen)  or  oxygen- 
deficient  (they  contain  molecular  hydrogen  or  carbon  monoxide 
or  both).  The  calculations  combine  a material  balance/equil- 
ibrium expression  such  as  the  one  obtained  by  combining  Equa- 
tions 6 and  11.  The  number  of  moles  of  water  vapor  present 
is  expressed  as 

H20  - ^S-02i/(*n  + o2*)  (18) 


where 

aH  * aH  “ ^ ~ OH 

i.e.,  the  number  of  hydrogen  atoms  available  for  the  water 
vapor  equilibrium.  Such  combined  expressions  have  mathematical 
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advantages  because  possible  negative  answers  are  avoided  for 
amounts  present. 

When  products  composition  has  been  determined,  a blast 
overpressure  can  be  computed  from  the  known  volume  of  the 
system  and  its  assigned  temperature  by  use  of  the  ideal  gas 
law.  This  overpressure,  however,  is  not  necessarily  the 
blast  maximum  peak  for  that  internal  explosion.  This  maxi- 
mum for  nonvented  explosions  corresponds  to  adiabatic  condi- 
tions where  the  energy  of  the  products  just  equals  that  for 
the  original  fuel  and  air.  Hence  to  ascertain  this  maximum, 
energy  considerations  become  involved. 


SOLID  CARBON  (SMOKE) 


Internal  explosions  that  are  very  fuel-rich  may  well 
show  nominal  products  with  elemental  carbon.  In  such  a cir- 
cumstance the  oxygen  present  is  not  sufficient  to  convert 
all  carbon  of  the  fuel  to  gaseous  products.  Particulate 
solid-phase  carbon  may  form  as  a smoke.  Systems  such  as 
these  with  a solid-phase  component  in  chemical  equilibrium 
with  gases  must  be  treated  in  a manner  that  differs  mathe- 
matically from  systems  of  gases  only  because  the  fugacity  of 
a solid  component  is  negligibly  affected  by  pressure.  Thus 
the  fugacity  does  not  appear  as  a variable  in  the  mathemati- 
cal expression  for  the  equilibrium  constant. 

The  pertinent  relations  for  the  smoke  reaction  are 
Formation  of  solid  carbon: 

2C0  =>  C02  + C Kn  - C02/C02  (19) 

This  mole  number  equilibrium  constant  Kn  can  be  evaluated 
from  the  thermodynamic  constants  of  formation  for  carbon  di- 
oxide and  for  carbon  monoxide,  perhaps  as  provided  by  the 
equations  of  Appendix  C.  A sample  calculation  is  shown  in 
Appendix  B. 

Calculations  for  the  possible  presence  of  solid  carbon 
in  equilibrium  with  gaseous  products  of  explosion  are  first 
computed  by  the  method  of  Appendix  B,  to  find  the  number  of 
moles  of  carbon  dioxide  and  carbon  monoxide  that  could  be 
present.  If  the  amount  of  carbon  dioxide  indicated  by  this 
computation  exceeds  that  prescribed  by  smoke  Equation  19 
(i.e.,  if  C02  >KnC0z),  no  solid  carbon  occurs  in  the  prod- 
ucts mixture.  However,  if  the  computed  number  of  moles  of 
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carbon  dioxide  is  less  than  that  prescribed  by  smoke  Equa- 
tion 19  (i.e.,  if  CO2  < ) > solid  carbon  may  be  present 
in  the  equilibrium  mixture.  In  this  circumstance  a trial 
value  for  the  number  of  moles  of  solid  carbon,  C,  is  selected. 
The  material  balance  for  elemental  carbon  (Equation  5)  is 
modified  to  the  form 

ac  - COo  + CO  + C (20) 

and  the  computation  illustrated  by  the  algorithm  of  Appendix 
D repeated.  Calculations  with  successive  trial  values  for 
solid  carbon  are  continued  until  conditions  of  smoke  Equation 
19  are  satisfied;  i.e.,  until  CO2  * KnC0z . This  second  order 
iteration  proceeds  readily  by  computer  where  successive  trial 
values  can  be  selected  automatically. 


CONSERVATION  OF  ENERGY  RELATIONS 


The  measure  for  energy  effects  in  a system  at  a speci- 
fied volume  is  the  system's  internal  energy,7  a thermodynamic 
item  that  differs  from  the  energy  item  usually  encountered — 
the  enthalpy  ("heat  content").  Two  rather  different  types 
of  internal  energy  are  involved;  one  thermal  in  nature,  the 
other  chemical.  The  thermal  aspect  is  associated  with  tem- 
perature and  heat  capacity;  the  chemical  aspect  is  associ- 
ated with  chemical  composition.  These  two  diverse  aspects 
can,  however,  be  combined  into  a single  term  if  they  relate 
to  a common  basis.  The  basis  selected  for  computations  is 
the  elements  in  their  ordinary  form  at  a reference  tempera- 
ture of  25®C  and  pressure  of  the  standard  atmosphere  (but 
pressure  level  is  relatively  immaterial  for  thermal  items). 

On  this  basis,  an  element  is  assigned  a chemical  energy  of 
zero,  and  its  energy  is  given  rather  simply  as  the  tempera- 
ture integral  of  its  heat  capacity.  A compound  such  as  car- 
bon dioxide  (or  an  element  not  in  its  ordinary  form  such  as 
monatomic  oxygen)  shows  both  this  thermal  aspect  and  an  ad- 
ditional chemical  aspect  known  as  energy  of  formation.  The 
total  energy  for  such  a material  becomes  the  sum  of  these 
two  separate  aspects. 


7 P.  J.  Kiefer,  G.  F.  Kinney,  and  M.  C.  Stuart.  Prin- 
ciples of  Engineering  Thermodynamics . New  York,  Wiley  and 
Sons,  1958. 
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CHEMICAL  ENERGY 


The  chemical  energy  aspect  for  a constant  volume  system 
.Ls  given  by  the  internal  energy  of  formation,  which  by  def- 
inition is  the  internal  energy  change  for  the  formation  re- 
action. This  item  is  not  directly  available  for  most  mate- 
rials. It  can,  however,  be  obtained  indirectly  by  standard 
thermochemical  methods  from  data  such  as  a heat  of  combus- 
tion as  provided  in  standard  sources  for  conventional  fuels 
and  for  explosives.3 

Several  items  with  regard  to  these  thermochemical  com- 
putations are  as  follows: 

1.  Thermochemical  data  commonly  relate  to  1 mole  of 
a material  with  definite  chemical  formula.  If  only  a chem- 
ical analysis  is  available,  it  is  convenient  to  express  this 
in  alternative  form  as  an  empirical  chemical  formula. 

2.  Many  published  therraochemical  data  follow  the 
older  algebraic  convention  where  heat  energy  leaving  a system, 
as  from  a combustion,  is  considered  to  be  algebraically  pos- 
itive. By  modern  thermodynamic  convention  when  the  energy 

of  the  system  is  diminished  by  such  an  exothermic  reaction, 
this  diminution  effect  is  algebraically  negative.  Hence  to 
utilize  published  thermochemical  data  for  thermodynamic  cal- 
culations, particularly  for  data  from  older  sources,  a change 
in  algebraic  sign  may  be  required. 

3.  Many  published  thermochemical  data  relate  directly 
to  combustion  reactions.  To  use  these  data  for  the  calcula- 
tions, they  first  must  be  converted  to  formation  values.  The 
calculation  is  a conventional  one  illustrated  in  Appendix  E. 
The  calculation  must  take  into  account  whether  water  formed 
in  the  combustion  is  considered  to  be  a vapor,  given  a 
"lower"  or  "net"  heating  value  (LEV),  or  a liquid,  giving  a 
"higher"  or  "gross"  heating  value  (HHV). 

4.  The  therraochemical  kilocalorie  is  used  as  an 
energy  unit  in  older  theraochemical  tables.  The  multiplying 
factor  for  conversion  into  joules  is  4184  exactly. 

5.  Most  tabulated  thermochemical  values  are  actually 
for  an  enthalpy  and  described  as  being  for  constant  pressure 
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conditions.  Conversion  of  these  iato  internal  energies  pro- 
ceeds readily  because  by  definition  the  enthalpy  is  the  sura 
of  The  internal  energy  and  the  pressure-volume  product  (Ap- 
pendix E).  Table  E-l  of  this  appendix  lists  internal  ener- 
gies of  formation  for  many  of  the  fuels  that  are  important 
here . 


6.  A formation  energy  represents  a small  difference 
between  large  numbers.  Its  numerical  value  is  thus  sensitive 
to  small  discrepancies  in  those  numbers.  These  small  dis- 
crepancies do  not  necessarily  contribute  substantial  uncer- 
tainty into  the  internal  blast  calculations. 

THERMAL  ENERGY 


The  thermal  component  of  an  internal  energy,  the  temper- 
ature integral  of  the  heat  capacity  at  constant  volume,  must 
relate  to  the  selected  reference  temperature  of  25°C.  The 
pertinent  interval  for  integration  then  becomes  that  from 
25aC  to  the  actual  temperature  for  fuel,  air,  or  explosion 
products,  respectively. 

Data  on  heat  capacity  for  many  fuels,  particularly  ones 
in  a condensed  phase,  may  not  be  readily  available.  Heat 
capacity  may  often  be  adequately  approximated  by  the  Kopp 
rule  for  additive  atomic  constants.  These  constants  have 
boen  evaluated  for  metric  units  and  are  given  in  Appendix  F. 
The  values  pertain  to  temperatures  near  room  temperature, 
and  within  their  inherent  uncertainties  apply  to  both  con- 
stant volume  and  const  \nt  pressure  conditions. 

Ordinary  air  contains  78-mol-%  nitrogen,  21-mol-%  oxygen , 
and  l-mol-%  argon.  It  is  a diatomic  ideal  gas  and  at  temper- 
atures near  258C  its  heat  capacity  at  constant  volume  has  the 
classic  value  of  Cv  a 5/2  R.  In  metric  units  the  gas  law  con- 
stant R is  8.31434  J/Craol^K),  so  that  the  molar  heat  capacity 
for  air  at  constant  volume  becomes 

Cv  - 5/2  x 8.314  “ 20.8  J/(mol-K)  (21) 

The  number  of  moles  of  air,  n , can  be  found  by  the  ideal  gas 
law  with  n * PV/RT.  Since  no  chemical  item  is  involved,  the 
total  internal  energy  for  air  at  temperature  T in  kelvins 
and  relative  to  the  elements  at  25"C  is 

2°  -20.8  (T  - 298 )n  (22) 

The  gaseous  products  from  combustion-explosion  reactions 
show  heat  capacities  that  are  conveniently  described  by  three 
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coefficient  equations.9  These  thermal  items  can  be  combined 
with  pertinent  chemical  items,  the  internal  energies  of  for- 
mation, to  give  a four-coefficient  equation  for  the  (total) 
internal  energy  for  each  component,  relative  to  the  elements 
at  25 °C . Thus 

£•*  * a + bT  + o?2  + d/T  (23) 

where  5*  is  the  standard  internal  energy  per  mole  at  temper- 
ature T and  a , b,  a,  and  d are  arbitrarily  selected  coeffi- 
cients. The  four  coefficients  for  this  equation  have  been 
evaluated  for  the  12  components  from  data  of  the  JANAF  tables 
for  the  temperature  range  500-4000  K.6  These  are  tabulated 
in  Appendix  G.  The  total  energy  for  a products  mixture  is 
then  a summation  of  individual  values  per  mole  times  the 
number  of  moles  of  each  component  present. 

ADIABATIC  COMBUSTIONS 


The  total  energy  of  the  products  (sum  of  chemical  plus 
thermal)  just  equals  that  for  the  original  fuel  plus  air  at 
the  adiabatic  flame  temperature  of  a combustion-explosion 
reaction.  A trial  value  for  products  temperature  is  assumed 
for  computation  of  this  flame  temperature.  Then  by  itera- 
tion, as  described  above,  an  associated  equilibrium  compo- 
sition and  internal  energy  of  products  for  this  assumed  tem- 
perature are  computed.  Successive  trial  values  for  tempera- 
ture are  selected  until  the  desired  equality  for  the  internal 
energies  is  found. 

Given  the  adiabatic  flame  temperature,  peak  products 
pressure  (for  a nonvented  internal  explosion)  can  be  computed 
by  the  ideal  gas  law  PV  » ( m/M)R7  where  m is  the  total  mass 
of  fuel  plus  air,  and  M the  formula  mass  for  the  products. 
Blast  overpressure  is  then  obtained  by  subtracting  the  am- 
bient atmospheric  pressure. 

INTERNAL  EXPLOSIONS  WITH  TNT 


The  conventional  explosive  TNT  (formula  C7H5N3O0)  also 
serves  as  a reference  for  explosives  in  general.  In  an 


9 C.  G.  Maier  and  K.  K.  Kelley.  "An  Equation  for  the 
Representation  of  High-Temperature  Heat  Capacity  Data,"  J. 
A.ner.  Chen.  3oa.  , Vol.  54  (1932),  p.  3243. 
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internal  explosion,  TNT  acts  as  an  explosive  fuel  where  its 
inherent  oxygen  content  supplements  that  of  the  air  and  so 
may  supply  a substantial  proportion  of  the  oxygen  for  the 
combustion-explosion  reaction.  Calculations  for  internal 
explosions,  as  outlined  above,  readily  permit  this;  the  re- 
sults of  such  calculations  are  reported  below. 

OVERPRESSURES 


Overpi .ooures  generated  in  internal  explosions  with  TNT 
were  shown  graphically  in  Figure  2 for  overpressures  in  bars 
(10s  pascals)  over  a three-decade  range  of  TNT-air  ratios  in 
kilograms  of  TNT  per  cubic  metre  of  air  at  25®C  and  1 bar. 
Also  shown  were  the  experimental  results.3  The  agreement 
obtained  is  gratifying  and  leads  to  confidence  in  the  thermo- 
dynamic computations. 

Both  calculated  and  experimental  results  for  internal 
explosions  with  TNT  show  a point  of  inflection  at  an  inter- 
mediate fuel-air  ratio,  although  this  was  not  recognized  in 
the  original  report  of  these  experiments.  Beyond  this  point 
fuel-air  ratios  are  such  that  the  primary  reaction  is  a det- 
onation where  increasing  overpressures  result  largely  from 
formation  of  additional  gases.  Below  this  inflection  point 
the  air  provides  sufficient  oxygen  so  that  combustion  effects 
in  the  afterburn  contribute  to  the  overall  temperature  rise 
and  resulting  pressure  increase. 

Four  stoichiometric  points  for  combustions  of  TNT  in  air 
are  shown  in  Figure  2.  The  overpressure  curve  proceeds 
smoothly  through  these  points.  The  CO-H2  point  corresponds 
to  the  most  fuel-rich  and  is  also  termed  the  smoke  point  be- 
cause with  a greater  amount  of  fuel  solid  carbon  may  be  one 
of  the  combustion  products. 

FLAME  TEMPERATURES 


Computations  for  internal  blast  overpressure  also  allow 
for  other  characteristics  of  the  blast.  Figure  3 shows  adia- 
batic flame  temperatures  for  the  combined  combustion-explo- 
sion reaction  of  TNT  in  air  at  25°C  and  1 bar  for  a wide 
range  of  TNT-air  ratios.  These  temperatures  are  appreciably 
higher  than  those  for  constant  pressure  conditions. 

In  contrast  with  pressures  for  internal  blast,  which 
increase  monotonically  with  fuel-air  ratio,  the  blast  tem- 
perature curve  for  TNT  reaches  a maximum  at  about  3055  K at 
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MASS-VOLUME  RATIO,  kg/m3 

FIGURE  3.  Temperatures  ta  Internal  Blast 
With  TNT-Air  Explosions. 


a TNT-air  ratio  of  about  0.7  kg/m3.  This  ratio  lies  between 
the  CO2-H2O  and  the  CO-H2O  stoichiometric  reference  points 
for  TNT  and  is  close  to  the  fuel-air  ratio  at  the  point  of 
inflection  in  the  overpressure  curve  This  maximum  corres- 
ponds to  maximum  utilization  of  oxygen  in  the  afterburn. 

PRODUCTS  FORMULA  MASS 

Chemical  dissociations  at  high  temperatures  give  pro- 
ducts such  as  monatomic  oxygen  and  hydrogen  and  reduce  the 
overall  formula  mass  for  the  products  mixture  (Figure  4). 
Formula  mass  versus  products  temperature  is  plotted  for  a 
representative  TNT  explosion  in  air — one  with  a fuel-air 
ratio  that  corresponds  approximately  to  maximum  flame  tem- 
perature. The  formula  mass  for  these  products  decreases  as 
much  as  10%  at  high  temperatures.  The  corresponding  increase 
in  the  number  of  moles  or  gases  present  can  make  an  appre- 
ciable contribution  to  blast  overpressures. 
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FIGURE  4.  Formula  Mass  for  Produ  ..s  From 
a TNT-Air  Internal  Explosion. 


PRODUCTS  COMPOSITION 

The  extent  of  chemical  dissociation  in  the  products  of 
an  internal  explosion  with  TNT  is  shown  in  Figure  5.  Com- 
position of  products  versus  the  TNT-air  ratio  (logarithmic 
scales)  is  plotted.  Maximum  dissociation  occurs  in  the  gen- 
eral region  of  the  CO2-H2O  and  the  CO2-H2  stoicioraetric  ref- 
erence points  and  at  approximately  the  fuel-air  ratio  for 
peak  adiabatic  temperature.  In  general  any  particular  pro- 
duct of  dissociation  is  present  only  in  a small  amount,  but 
together  these  dissociations  have  substantial  influence  on 
products  temperature  because  of  their  pronounced  endothermic 
nature. 

PRODUCTS  HEAT  CAPACITY  AND  HEAT  CAPACITY  RATIO 

Heat  capacity  data  are  important  for  many  explosion  cal- 
culations. The  heat  capacity  for  a system  at  constant  volume 
is  the  temperature  derivative  of  its  internal  energy,  which 
can  be  found  directly  from  results  of  the  energy  calcula- 
tions outlined  above.  The  molar  heat  capacity  Cv  (in  joules 
per  mole  kelvin)  has  been  computed  for  the  products  of  the 
explosion  of  1 kg  of  TNT/m3  of  air  at  25“C  and  1 bar.  These 
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TNT-AIR  RATIO,  kg/m3 

FIGURE  5.  Composition  of  Products  for  Internal 
Explosions  With  TNT-Air  Mixtures. 

lata  are  presented  as  a function  of  products  temperature  in 
Figure  6.  This  heat  capacity  is  greatly  influenced  by  chem- 
ical dissociation  effects  and  increases  monotonically  with 
temperature.  At  higher  temperatures  it  appreciably  exceeds 
both  the  conventional  diatomic  value  of  3/2 R =»  20.3  J/raol-K) 
and  the  triatomic  value  of  7/2 R * 29.1  J/(mol-K). 

Calculations  such  as  those  for  compressible  fluid  flow 
may  utilize  heat  capacity  data  in  the  form  of  a ratio  of 
constant  pressure  to  constant  volume.  For  ordinary  gases  a 
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TEMJERATUR*.  K 


FIGURE  6.  Molar  Heat  Capacity  Cv 
for  Products  From  an  Internal  TNT- 
Air  Explosion. 


simple  relationship  may  not  apply  here  because  of  the  effects 
of  chemical  dissociation  on  formula  mass.  A corresponding 
ratio  k » (Cv+fl)/Cy  ersus  temperature  is  plotted,  however, 
in  Figure  7 for  the  products  for  three  different  TNT-air 
ratios.  The  ratio  decreases  monotonically  with  temperature 
in  each  instance,  but  a changing  products  composit ion  appears 
o prevent  any  simple  correlation  between  heat  capacity  ratio 
and  fuel-air  ratio. 

OTHER  EXPLOSIVE  FUELS 

Explosive  fuels  generally  show  internal  explosion  over- 
pressures similar  to  those  for  TNT — overpressures  that  in- 
crease monotonically  with  fuel-air  ratio.  The  more  highly 
oxygenated  fuels  require  less  oxygen  from  the  air  for  com- 
plete reaction  and  consequently  show  lesser  afterburn  effects. 
Figure  8 shows  overpressures  for  reference  internal  explo- 
sions with  TNT,  an  oxygen-deficient  explosive  (oxygen  balance 
of  -74%  to  CC>2  and  -25%  to  CO)8  with  pronounced  afterburn 
effects  that  depend  on  the  TNT-air  ratio.  Shown  for  compari- 
son are  the  overpressures  for  internal  explosions  with 
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FIGURE  7.  Heat  Capacity  Ratio  k * 

( Cv  + R)/C,j  for  Explosion  Products. 


pentaerythritol  tetranitrate  (PETN) , an  explosive  approxi- 
mately balanced  in  oxygen  (oxygen  content  between  the  CO^- 
H2O  and  the  CO-HgO  stoichiometric  reference  points  and  with 
oxygen  balances  of  -10%  to  CO2  and  +15%  to  CO).  Because  a 
small  afterburn  effect  is  possible,  only  a small  inflection 
point  is  indicated  in  the  overpressure  curve.  Also  shown  is 
the  overpressure  curve  for  the  oxygen-rich  explosive  mannitol 
hexanitrate  (MHN ) with  more  oxygen  even  than  that  for  the 
CO2-H2O  stoichiometric  reference  point,  and  with  an  oxygen 
balance  of  +7%  to  CO2  and  +28%  to  CO.  No  afterburn  is  pos- 
sible, and  there  is  no  point  of  inflection  in  the  overpres- 
sure curve. 

Overpressures  developed  in  internal  explosions  with 
many  explosive  fuels,  both  oxygen-rich  and  oxygen-deficient, 
have  been  computed  for  a wide  range  of  fuel-air  ratios  with 
air  initially  at  1 bar  and  25’C  (Tables  1 and  2).  For  these 
explosive  fuels  the  overpressure  always  increases  with  the 
fuel-air  ratio,  but  actual  overpressures  depend  on  the  type 
of  fuel.  A characterization  of  this  aspect  of  internal 
blast  follows. 

INTERNAL  BLAST  YIELDS 

In  the  low  overpressure  range,  the  overpressures  devel- 
oped by  a given  fuel-air  ratio  of  PETN  are  less  than  those 
for  TNT.  The  internal  blast  yield  for  PETN  is  therefore 
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FUEL-AIR  RATIO,  kg/m3 

FIGURE  8.  Overpressures  for  Internal  Explosions  With 
Three  Explosives  Having  Different  Oxygen  Content. 

defined  to  be  less  than  that  for  TNT  (Figure  8).  Such  in- 
ternal blast  yields  are  conveniently  expressed  relative  to 
TNT,  an  explosive  with  a definite  chemical  'ormula  and  re- 
producible characteristics.  The  internal  blast  yield  for 
PETN  in  the  low  overpressure  range  is  only  52%  of  that  for 
TNT  on  a mass  basis;  i.e.,  52  grams  of  TNT  give  the  same  in- 
ternal blast  overpressure  as  100  grams  of  PETN  (in  the  low 
overpressure  range).  On  a volume  basis  where  relative  den- 
sity must  be  taken  into  account,  the  yield  for  PETN  is  56% 
of  that  for  TNT'  i.e.,  56  ctn3  of  TNT  give  the  same  internal 
blast  as  100  cm3  of  PETN  (in  the  low  overpressure  range, 
which  corresponds  to  explosions  in  a relatively  large  volume 
of  air). 
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Internal  blast  yields  can  readily  be  assigned  by  a 
graphical  method  that  utilizes  plots  of  blast  overpressure 
versus  logarithm  of  the  fuel-air  ratio.  The  plots  for  the 
two  materials  are  superimposed  so  that  overpressures  in  the 
range  of  interest  closely  coincide.  The  relative  internal 
blast  yield  for  this  range  of  overpressures  is  then  given  as 
the  inverse  ratio  of  the  two  corresponding  fuel-air  ratios. 

Blast  yields  relative  to  TNT  for  a large  number  of  ex- 
plosive fuels  are  given  in  Table  3.  Two  values  are  included 
for  each  of  these  explosive  fuels,  one  for  the  low  overpres- 
sure range  (on  the  order  of  5 bars)  and  one  for  high  over- 
pressures (20  to  30  bars).  These  two  values  may  be  quite 
different . 

In  the  low  overpressure  range,  the  yield  includes  con- 
tributions from  the  afterburn  effect,  which  can  be  quite 
significant  for  oxygen-deficient  explosives.  At  low  over- 
pressures an  oxygen-deficient  explosive  fuel  with  substan- 
tial afterburn  effects  (dinitrotoluene , for  example),  can 
generate  a more  damaging  internal  explosion  than  an  oxygen- 
rich  explosive  such  as  mannitol  hexanitrate. 

The  higher  overpressures  generated  in  internal  explo- 
sions with  explosive  fuels  result  primarily  from  detonations, 
and  so  are  little  influenced  by  the  afterburn.  An  explosive 
fuel  with  high  oxygen  content  does  not  utilize  all  its  oxy- 
gen, and  the  excess  oxygen  is  simply  a diluent  that  reduces 
the  yield  in  the  higher  overpressure  range.  Maximum  internal 
blast  yields  in  the  high  overpressure  range  are  obtained  from 
explosive  fuels  that  are  approximately  balanced  in  oi.ygen. 

The  pronounced  effect  of  the  oxygen  content  of  an  explo- 
sive on  its  internal  blast  yield  is  illustrated  in  figure  9. 
The  overpressure  in  both  the  low  and  high  ranges  are  shown 
for  the  seven  members  of  the  benzene  series  of  nitrated  ma- 
terials. Neither  benzene  (CgHg)  nor  mononitrobenzene 
(C6H5NO2)  detonates,  and  hence  they  do  not  achieve  the  high 
overpressure  range.  They  do,  however,  show  appreciable 
overpressures  in  the  low  overpressure  range  as  generated  by 
combustion  reactions.  In  this  range  these  particular  over- 
pressures are  substantially  higher  than  those  for  the  highly 
nitrated  materials  with  their  greater  proportion  of  noncom- 
bustible constituents. 

In  the  high  overpressure  range,  materials  with  either 
low  oxygen  or  high  oxygen  content  lead  to  lower  overpressures 
than  the  overpressures  obtained  from  the  same  amount  of  an 
oxygen-balanced  material.  To  illustrate  for  the  seven-member 
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FIGURE  9.  Overpressures  for  the 
Nitrated  3enzene  Series. 
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nitrated  benzene  series,  tetraaitrobenzena  (an  explosive 
approximately  balanced  for  the  CO-H2O  stoichiometric  refer- 
ence point)  is  most  effective  from  a yield  viewpoint  at  high 
overpressure  (Figure  10). 


INTERNAL  EXPLOSIONS  WITH  CONVENTIONAL  COMBUSTIBLES 


Overpressures  generated  in  the  internal  blast  for  three 
conventional  fuels — benzene,  JP-4  (a  hydrocarbon  fuel  with 
an  empirical  formula  C9H17),  and  ethylene  oxide — are  plotted 
as  a function  of  fuel-air  ratio  (Figure  11).  The  correspond- 
ing curve  for  the  explosive  fuel  TNT  is  included  for  compari- 
son. All  the  curves  for  three  conventional  combustibles  show 
a maximum  overpressure  on  the  order  of  9 to  11  burs  that 
characteristically  lies  between  the  CO2-H2O  and  the  CO-H2O 
stoichiometric  reference  points.  The  two  hydrocarbon  fuels, 
benzene  and  JP-4,  show  typical  behavior  even  though  their 
hydrogen-carbon  ratios  are  quite  different.  Ethylene  oxide 
(C2H2O) , an  oxygen-containing  fuel,  also  shows  a maximum  in 
its  overpressure  curve,  but  one  that  is  broader  and  that 
occurs  at  an  appreciably  greater  fuel-air  ratio.  This  is  a 
consequence  of  the  diluting  effect  of  oxygen  in  the  molecule. 
For  some  purposes  the  greater  mass  of  ethylene  oxide  required 
to  achieve  a maximum  blast  overpressure  can  be  a disadvantage, 
while  the  broad  range  for  high  overpressures  can  be  an  ad- 
vantage. 

The  internal  blast  behavior  of  selected  conventional 
fuels  is  summarized  in  Tables  4 and  5 for  initial  conditions 
with  air  at  1 bar  and  25"C.  The  maximum  overpressure  that 
each  fuel  can  generate  in  confined  air  is  shown  along  with 
the  required  fuel-air  ratio.  The  tables  also  show  maximum 
flame  temperatures  and  the  associated  fuel-air  ratio.  In 
general  the  fuel-air  ratio  for  maximum  temperature  is  approx- 
imately that  for  maximum  blast  overpressure.  This  is  to  be 
expected  because  at  lesser  ratios  some  of  the  available  oxy- 
gen is  not  utilized,  while  at  greater  ratios  any  excess  fuel 
contributes  only  to  the  mass  of  products  to  be  warmed  and 
gives  reduced  temperature  and  pressure  rises. 

The  flame  temperatures  given  in  Tables  4 and  5 for  com- 
bustions at  constant  volume  are  appreciably  greater  than 
those  for  constant  pressure  conditions:  (1)  the  beat  capa- 
city at  constant  volume  is  less  than  that  at  constant  pres- 
sure so  that  a given  energy  release  produces  a greater  tem- 
perature rise,  and  (2)  the  chemical  dissociations  in  the 
flame  are  partially  suppressed  by  the  higher  pressure  condi- 
tions of  constant  volume. 
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NUMBER  OF  NITRO  GROUPS  PER  BENZENE  NUCLEUS 


FIGURE  10.  Internal  Blast  Yields 
for  the  Nitrated  Benzene  Series. 


GASEOUS  FUELS 


A confined  mixture  of  a gaseous  fuel  such  as  methane  in 
air  can  explode  only  when  its  fuel-air  ratio  is  within  cer- 
tain limits.10  For  a mixture  of  methane  in  air,  the  lower 
flammability  limit  is  5 tnol%  (5  vol%)  and  the  vapor  flam- 
mability limit  is  15%.  Only  within  these  limits  can  a mix- 
ture of  methane  in  air  produce  an  internal  explosion.  Adia- 
batic values  for  the  internal  blast  overpressures  computed 
for  these  limiting  ratios  are  5.5  and  7.0  bars,  respectively. 
Maximum  overpressure  for  a methane-air  explosion  is  8.0  bars, 
which  occurs  at  10.0  vol%  methane,  a composition  intermedi- 
ate between  the  flammability  limits,  and  also  between  those 


10  Bernard  Lewis  and  G.  von  Elbe.  Combustion,  Flames, 
and  Explosions  of  Gases.  New  Yorh,  Academic  Press,  1951. 
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FIGURE  11.  Blast  Overpressures  Versus  Fuel-Air 
Ratio  for  Three  Nonexplosive  Fuels. 


for  the  CO2-H2O  (9.5%)  and  the  CO-H2O  (12.3%)  stoichiometric 
points.  Calculations  for  the  internal  blast  from  various 
gaseous  fuels  are  based  on  initial  conditions  of  25®C  and  1 
bar  total  pressure  for  the  mixture  (Table  5). 

Table  5 also  provides  the  corresponding  internal  blast 
yields  expressed  in  terms  of  methane  for  gaseous  fuels.  The 
yield  is  the  volume  of  a gas  (moles)  or  the  mass,  relative 
to  methane,  required  to  produce  the  same  overpressure  upon 
its  explosion  in  air.  Such  yields  vary  with  the  nature  of 
the  gas,  with  values  particularly  sensitive  to  oxygen  con- 
tent of  the  gas  molecule.  Yields  in  terms  of  volume  are 
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appreciably  greater  for  the  more  dense  gases  than  for  the 
lighter  gases. 

FLAME  EMISSIVITY 

Overpressures  reported  in  Tables  4 and  5 were  computed 
for  adiabatic  conditions;  i.e.,  they  do  not  provide  for  heat 
transfer  from  hot  explosive  products  to  cool  confining  walls 
either  by  radiant  energy  or  by  conduction-convection  pro- 
cesses. Heat  transfer  effects  are  very  time-dependent  and 
so  are  more  important  for  the  relatively  slow  explosions 
with  conventional  fuels  than  for  the  rapid  ones  with  explo- 
sive fuels.  With  an  explosive  fuel  the  detonation  is  very 
rapid  indeed  and  also  serves  to  accelerate  combustion  ef fects 
in  the  afterburn  because  of  an  increased  temperature.  The 
relatively  minor  effect  of  heat  transfer  with  explosive  fuels 
is  indicated  in  Fig  re  2,  where  the  peak  overpressures  com- 
puted for  TNT  on  an  adiabatic  basis  agree  quite  well  with 
those  observed  experimentally. 

The  importance  of  heat  transfer  iu  the  relatively  slow 
internal  explosions  with  conventional  combustibles  exploding 
in  coal  mines  or  in  grain  elevators  can  be  deduced  from  ex- 
perimental data,11  which  data  are  compared  in  Figure  12  with 
results  of  the  adiabatic  calculations  made  as  described 
above.  Figure  12(a)  is  for  coal  dust  in  air,  Figure  12(b)  for 
starch  in  air;  and  Figure  12(c)  for  methane  in  air  for  two 
conditions,  quiescent  mixtures  and  initially  turbulent  mix- 
tures. Figure  12(c)  illustrates  the  effp-  of  rate  of  com- 
bustion on  peak  overpressure  developed  in  : ternal  explo- 
sions. Turbulent  mixtures  with  their  morr‘  rapid  rata  of 
combustion  have  less  time  for  heat  transfer,  and  there  is  a 
corresponding  lesser  reduction  in  peak  overpressure  from  the 
adiabatic  values.  The  cooling  effect  of  confining  walls  can 
nevertheless  cause  a 10  tc  2C%  reduction  in  peak  blast  over- 
pressure from  explosions  with  conventional  combustibles. 

These  effects  of  heat  transfer  can  be  expressed  quanti- 
tatively in  terras  of  the  amount  of  energy  transferred  during 


llU.S.  Bureau  of  Mines.  explosion  Development  in 
Closed  Vessels,  by  John  Nagy,  E.  C.  Seiler,  J.  W.  Conn,  and 
H.  C.  Varakis.  Washington,  D.C.,  U.S.  Bur.  Mines,  1971. 
(Bureau  of  Mines  Report  of  Investigations  7507,  publication 
UNCLASSIFIED. ) 


fflic  TP  6089 


the  explosion  process,  This  computation  finds,  by  iterative 
methods,  the  products  temperature  and  its  associated  internal 
energy  that  correspond  to  the  observed  overpressure.  These 
are  then  compared  with  adiabatic  values  as  computed  by  meth- 
ods described  above. 

The  important  mechanism  for  heat  transfer  at  high  tem- 
perature is  thermal  radiation,  and  its  associated  coefficient 
is  the  overall  flame  emissivity.  This  overall  emissivity  is 
a tenuous  item  that  depends  on  factors  such  as  flame  temper- 
ature, spectral  distribution  of  the  energies  present,  path 
length  for  the  radiant  energy,  and  on  possible  conduction- 
convection  effects.  An  overall  emissivity  can  still  be  as- 
signed in  any  particular  circumstance  when  the  rate  of  heat 
flow,  plus  the  emitting  and  absorbing  areas  and  their  temper- 
atures, is  known.  Such  is  the  case  for  several  of  the  ex- 
perimental explosions  reported.11  Appendix  H illustrates  a 
rep-esentat ive  computation. 

The  results  of  a series  of  emissivity  calculations  are 
given  in  Table  6.  Appreciable  scatter  is  observed  and  re- 
flects experimental  difficulty  in  providing  reproducible 
conditions.  Tentative  conclusions,  however,  can  be  deduced. 
Turbulent  flames  show  appreciably  higher  overall  eraissivities 
than  do  quiescent  flames,  with  average  emissivity  values  of 
0.55  and  0.22,  respectively.  This  enhanced  emissivity  of  the 
turbulent  flame  can  be  correlated  with  the  observation  that 
turbulent  flames  are  more  luminous  and  substantially  thicker 
than  quiescent  ones  and  are  therefore  more  effective  thermal 
emitters.  The  table  also  indicates  that  confining  vessel 
volume  has  a lesser  effect  than  does  the  degree  of  turbulence. 


INJECTION  VELOCITY  EFFECT 


The  kinetic  energy  of  injection  must  be  added  to  the 
system  when  fuel  for  an  internal  explosion  is  injected  into 
a confining  volume.  Figure  13  shows  the  effect  of  injection 
energy  on  cyclohexane  explosions.  A gain  in  overpressure  of 
0.7  bar  is  realized  for  a fuel-air  ratio  of  0.1  kg/m3  (opti- 
mum for  sea  level)  as  the  velocity  of  fuel  injection  is  in- 
creased from  0 to  the  rather  substantial  3000  m/s,  a veloc- 
ity that  corresponds  to  a Mach  number  of  8.7  in  ordinary  air. 
Thus  it  appears  that  the  injection  effect  is  not  important 


except  for  large  injection  velocities.  The  system  is  over- 
fueled at  a greater  fuel-air  ratio  (0.15  kg/m’),  but  an  in- 
crease in  the  injection  velocity  from  0 to  3000  m/s  gives  a 
larger  overall  increase  in  overpressure  cf  1.8  bars. 


38 


ymi  mull,  mum  uiiiiiii  miiiiiiiii  uiimLiC*...  .iWUwUi.ai  to..  iUHU[..,mil  y ...m  i^lllOfc.HUnutLiiiiilll  IlllllMlI  llfaidlilllhlllllly  M lUiiilidli  Mum  luiL.'tilljHJJUU  UllUllIlilll  ilil^ 


nidtiljliiiliilliULlllJn  inlM'lllliiUilifcJ  Jl 


NNC  TP  6089 


FIGURE  13.  Increased  Velocity  of  Fuel 
Injection  Causes  Moderate  Increase  in 
Blast  Overpressure. 


ALTITUDE  EFFECT 

An  increase  in  altitude  reduces  the  amount  of  oxygen  in 
the  system  when  a fuel  burns  in  a confined  volume  of  air  at 
ambient  pressure.12  For  a constant  amount  of  fuel,  this 


1 2 Handbook  of  Tables  for  Applied  Engineering  Soienae , 
2nd  ed.  "U.S.  Standard  Atmosphere,  1962,"  Cleveland,  Ohio, 
CRC  Press,  1976.  P.  651. 
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causes  a shift  toward  richer  mixtures  and  a reduction  in  the 
maximum  overpressure  that  can  be  developed.  This  decrease 
in  maximum  overpressure  is  shown  in  Figure  14  for  TNT  and  a 
conventional  hydrocarbon  fuel,  with  both  at  a fuel-air  ratio 
of  92  g/o3,  a ratio  that  is  near  the  optimum  for  the  hydro- 
carbon fuel. 


The  effect  of  the  reduction  in  overpressure  with  alti- 
tude can  be  somewhat  compensated  for  by  varying  the  amount 
of  fuel  injected  to  keep  the  fuel-air  ratio  near  the  optimum 
(Figure  14).  This,  however,  represents  a lesser  total  amount 
of  fuel  and  air  present  so  that  again  the  overpressure  de- 
creases sharply  with  altitude. 


FIGURE  14.  Inte-^al  Blast  Overpressure 
Decrease  at  Inc. eased  Altitude. 
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SUMMARY 


Confined  explosions  sued  as  those  in  grain  elevators, 
coal  mines,  or  ship  compartments  can  be  very  damaging.  The 
peak  overpressures  that  such  explosions  can  generate  have 
been  computed  from  first  principles  of  conservation  of  mass, 
conservation  of  energy,  and  minimum  Gibbs  free  energy.  The 
amounts  of  the  various  chemical  species  that  are  present, 
including  those  resulting  from  chemical  dissociation  in  the 
products  mixture  at  explosion  temperature,  are  found  by  using 
iterative  methods.  This  iteration  is  readily  performed  on  a 
small  tabletop  computer.  The  peak  overpressure  is  then  found 
by  the  ideal  gas  law  from  products  composition  and  tempera- 
ture. Other  items  such  as  formula  mass  and  heat  capacity  for 
the  products  can  also  be  determined.  Computer  overpressures 
agree  closely  with  experimental  observations  on  confined  ex- 
plosions of  TNT  in  air  for  a three-decade  range  of  fuel-air 
ratios. 

Peak  overpressures  have  been  computed  for  a large  number 
of  confined  explosions  with  explosive  fuels,  both  pure  chem- 
ical explosives  and  mixtures  of  these,  and  for  a wide  range 
of  fuel-air  ratios.  They  pertain  directly  to  nonvented  adia- 
batic explosions  that  occur  so  rapidly  that  there  is  little 
time  for  gas  leakage  or  for  heat  transfer  effects.  These 
overpressures  increase  monotonically  from  very  low  values  at 
low  fuel-air  ratios  up  to  more  than  70  bars  at  high  fuel-air 
ratios  (Tables  1 and  2). 

The  relative  effectiveness  of  an  explosive  fuel  in  gen- 
erating blast  overpressure  can  be  described  in  terms  of  an 
internal  blast  yield.  This  is  the  inverse  ratio  of  the 
amount  of  an  explosive  relative  to  some  standard,  here  taken 
as  TNT,  that  is  required  for  producing  the  same  overpressure. 
These  yields  depend  on  both  tne  type  of  explosive  and  the 
fuel-air  ratio.  Thus  oxygen-deficient  explosives  in  the  low 
fuel-air  range  show  pronounced  afterburn  effects  that  con- 
tribute to  the  pressure  rise.  But  in  the  high  fuel-air  range, 
relatively  less  oxygen  is  available  and  afterburn  is  less  im- 
portant. Also,  for  oxygen-rich  explosives  afterburn  effects 
are  not  pertinent.  Considerations  such  as  these  make  it  pos- 
sible to  correlate  the  chemical  structure  of  an  explosive 
with  the  internal  blasts  it  can  generate.  These  considera- 
tions also  indicate  that  two  different  internal  blast  yield 
values  should  be  assigned  to  a given  explosive — one  yield 
value  for  the  low  overpressure  range  and  one  for  the  high. 
These  yield  pairs  are  given  in  Table  3 for  representative 
chemical  explosives  and  explosive  mixtures. 
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Internal  explosions  with  ordinary  combustibles  such  as 
fuel  oil,  coal  dust,  or  starch  can  also  be  important.  For 
these  explosions  the  maximum  adiabatic  overpressures  attained 
are  usually  on  the  order  of  8 bars,  at  a definite  fuel-air 
ratio  that  characteristically  lies  between  the  C02-H20  and 
the  CO-H2O  stoichiometric  reference  points.  At  fuel-air 
ratios  below  the  maximum  overpressure  there  is  unutilized 
oxygen  from  the  air,  and  above  it  the  extra  fuel  acts  only 
as  material  to  be  warmed,  depressing  the  temperature  rise  and 
associated  pressure  rise.  Maximum  adiabatic  overpressures, 
along  with  associated  flame  temperatures  and  fuel-air  ratios, 
are  given  for  a number  of  liquid  and  solid  combustibles  in 
Table  4 and  for  combustible  gases  in  Table  5.  Values  at  the 
lower  and  upper  explosive  limits  are  also  included  in  Table 
5.  Blast  yields  relative  to  benzene  (for  liquids  and  solids) 
or  relative  to  methane  (for  gases)  are  also  tabulated. 

All  peak  overpressures  reported  above  are  for  adiabatic 
explosions;  i.e.,  ones  without  heat  transfer  from  explosion 
flame  to  confining  walls.  They  pertain  directly  to  explosive 
fuels  where  the  detonation  and  afterburn,  if  any,  are  very 
rapid  indeed.  In  internal  explosions  with  conventional  com- 
bustibles however,  the  combustion  reaction  may  be  suffi- 
ciently slow  so  that  noticeable  heat  transfer  can  occur. 

Thus  for  methane  explosions  in  air,  observed  overpressures 
are  only  about  80%  of  those  computed  for  adiabatic  conditions. 
From  such  data  it  is  calculated  that  the  overall  thermal 
emissivity  for  quiescent  methane  flames  is  about  0.22,  and 
about  0.55  for  the  more  luminous  flames  from  a turbulent 
methane-air  mixture. 

The  calculated  overpressure  effects  of  fuel  injection 
at  high  velocity  are  quite  modest.  In  addition,  calcula- 
tions for  the  effect  of  altitude  on  internal  blast  show  a 
negative  trend,  i.e.,  lower  air  pressures  at  higher  alti- 
tudes reduce  the  amount  of  oxygeu  available  for  the  explo- 
sion reaction,  resulting  in  a substantial  decrease  in  over- 
pressure. 
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Appendix  A 

STOICHIOMETRIC  POINT  CALCULATIONS 


1.  Find  nominal  products  from  the  reaction  of  TNT  with 
air  to  the  CO-H2O  stoichiometric  reference  point  in  air  at 
25*C  and  1 bar. 


I 


[■ 


2.  What  is  the  corresponding  products  composition  in 
mole  percent,  and  the  corresponding  fuel-air  ratio  in  kilo- 
grams per  cubic  metre  (stere)?  The  formula,  for  TNT  is 
C785N3O0,  anc*  its  formula  ma3s  is  227.13  g/mol.  Air  is  21% 
oxygen,  78%  nitrogen,  and  1%  argon. 

Ansuera . Selecting  1 mole  of  TNT  as  a basis  and  ex- 
pressing  the  pertinent  combustion  reaction  in  the  form  of  a 
chemical  equation, 

C7H5N30b  + 1.75  02  - 7 CO  + 2.50  H20  + 1-50  Ng 

The  number  of  moles  of  air  corresponding  to  1.75  moles 
of  oxygen  is  1.75/0.21  » 8.33.  Nitrogen  supplied  by  this 
air  is  thus  8.33  x 0.78  - 6.50  moles,  to  which  is  added  1.50 
moles  from  the  TNT  reaction  to  give  a total  of  8.00  moles 
nitrogen  in  the  products.  Argon  is  8.33  x 0.01  =*  0.08  mole. 
Expressed  in  tabular  form, 


Component 

Moles 

Mole  % 

CO 

7.00 

39.8 

h2o 

2.50 

14.2 

n2 

8.00 

45 . 5 

Ar 

0.  08 

0.5 

Total 

17.58 

100.0 

The  volume  of  the  8.33  moles  of  air  in  this  reaction,  at  25°C  I 
and  1 bar,  is  given  by  the  ideal  gas  law  as  8.33  x 0.083  14  x | 
298.15  * 206.6  litres,  or  0.2066  m3.  The  mass  of  1 mole  of  j 
TNT  is  0.227  13  kilogram  so  that  the  fuel-air  ratio  becomes  1 
0.227  13/0.20 66  = 1.100  kg/m3.  Corresponding  calculations  for  ■] 
the  CO2-H0O  point  show  a fuel-air  ratio  of  0.367  kg/o3 , and  1 
for  the  CO2-H2  and  CO-H2  points,  ratios  of  0.481  and  3.848  a 
kg/m3,  respectively.  i 
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Appendix  B 

FREE  ENERGY  COMPUTATIONS 


1.  What  is  the  thermodynamic  equilibrium  constant  K 
for  the  dissociation  of  1 mole  of  carbon  dioxide  at  2500 
as  given  by  data  from  the  JANAF  tables?6 


Ans'Jar.  Expressing  the  dissociation  as  a chemical  equa- 
tion , 

C02  - CO  + 1/2  02 

its  thermodynamic  equilibrium  constant  Up  is 

(CO)  x (02)* 

xp  “ (cop 

where  parentheses  denote  the  fugacity  of  an  indicated  mate- 
rial (the  thermodynamic  item  fugacity  corresponds  to  a "par- 
tial pressure"). 

The  tables  provide  directly  the  (logarithm  to  base  ten 
of  the)  equilibrium  constant  of  formation  for  each  component 
of  the  equilibrium.  Such  values  are  to  be  combined  just  as 
are  the  fugacities  in  the  equilibrium  constant  expression. 

For  carbon  dioxide  at  2500  K,  the  formation  logarithm  is 
8.280  and  the  constant  itself  is  1.905  x 10s.  Corresponding 
values  for  carbon  monoxide  are  6.840  and  6.918  x 106 . For 
molecular  oxygen,  an  element  in  its  ordinary  form,  the  values 
are  zero  and  unity,  respectively.  Combining, 


6.918  x 106  x (1)4 
p 1.905  x 103 


* 3.63  x 10"2 


The  numerical  value  K*  - 3.63  x 10“  applies  at  the  specified 
temperature  of  2500  IT  with  fugacities  relative  to  the  arbi- 
trary standard  of  1.013  25  bars  (the  standard  atmosphere). 


2.  What  is  the  corresponding  mole  number  equilibrium 
constant  Kn  at  2500  K for  the  dissociation  of  carbon  dioxide 
in  a system  with  a volume  of  1 ra3? 
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4nawaj».  The  calculation  proceeds  by  noting  that  a fu- 
gacity  in  an  ideal  gas  mixture  is  given  as  the  product  of  a 
mole  fraction  and  system  pressure.  Thus  the  fugacity  for 
carbon  dioxide  (CO2)  * CO2  x P/n,  where  CO2  is  the  mole  num- 
ber for  the  carbon  dioxide  present,  P the  pressure  in  the 
system,  and  n the  total  number  of  moles  of  gases.  For  ideal 
gases  PV  ■ nRT,  and  the  ratio  P/n  is  also  given  as  RT/V  where 
R is  the  gas  law  constant,  T the  absolute  temperature,  and  V 
the  volume  of  the  system.  Substituting  into  the  expression 
for  the  thermodynamic  equilibrium  constant  and  rearranging, 


(RTJV)' 


3.63  x 10~2 


where  again  ?.  chemical  formula  indicates  the  number  of  moles 
of  the  specified  species  in  the  equilibrium  mixture.  The 
central  term,  with  mole  numbers  only,  also  defines  a mole 
number  equilibrium  constant  Xn , or 


CO  x Oo* 
— C0^ 


Kv(V/RT)b  « 3.63  x 102  (VfRT)k 


In  this  expression  the  item  (V/RT)  has  units  of  a fugacity 
at  the  standard  atmosphere.  Then  selecting  litres  for  volume, 
the  appropriate  value  for  R is  0.082  06  litre  ( 1 ■ atm) / (mol • K) , 
so  that  at  2500  K 

* 3.63  x 10"2  x (1000/0.082  06  x 2500)* 

8.017  x 10-2 

Such  mole  number  constants  greatly  facilitate  equilibrium 
calculations  for  systems  with  specified  volume  and  tempera- 
ture. 

3.  Compute  the  mole  number  equilibrium  constant  for  the 
smoke  reaction  in  a system  with  a volume  of  1 m3  at  2500  K. 

Answer ■ The  smoke  reaction  equilibrium  can  be  repre- 
sented as 


2 CO  * CO2  + C(S) 

Noting  that  the  fugacity  of  solid  carbon  is  unity,  the  ther- 
modynamic equilibrium  constant  becomes 


Xp  - (C02)/(C0): 


mump* 1 
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Thermodynamic  equilibrium  constants  of  formation  for  carbon 
dioxide  and  carbon  monoxide  at  2500  K are  Indicated  above  as 
1.905  x 10®  and  6.918  x 10s.  The  thermodynamic  constant  for 
the  smoke  reaction  becomes 

Xp  m (C02)/(C0)2  - 1.905  x 10®/(6. 918  x 106)2  - 
3.980  x 10"® 

This  thermodynamic  constant  Kp  is  converted  to  the  mole  num- 
ber constant  Kn  by  the  raethodillustrated  in  paragraph  2 
above , or 

Hn  - C02/C02  - Kp(RT/Y)  - Kp  X 0.082  06  x 2500/1000 
Kn  - 8.166  x 10W 
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Appendix  C 

EQUILIBRIUM  CONSTANTS  OF  FORMATION 


Logarithm  to  base  ten  of  the  ♦‘hermodynamic  equilibrium 
constant  of  formation  Kf  for  self  ;te d combustion  species  is 
presented  as  a function  of  absolute  temperature  T in  kelvins 
Values  are  based  on  ideal  gases  and  the  standard  state  of  1 
standard  atmosphere  ( 1 . 013  25  bars  ) ana  apply  over  the  tem- 
perature range  1500-4000  K. 

logr10  Kf  - 4 + S/(T  - C) 


TABLE  C— 1 . Coefficients  for  Calculating 
Equilibrium  Constants  of  Formation. 


Chemical 

Formula 

A 

B 

c 

co2 

-0.099 

21  096 

18 

CO 

4.069 

7 437 

184 

h2o 

-3.100 

13  415 

20 

0 

3 . 545 

-13  553 

16 

H 

3.242 

-12  277 

35 

OH 

0.662 

-1  723 

-112 

N 

3.582 

-25  435 

14 

NO 

0.645 

-4  652 

-15 
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Appendix  D 
THE  ALGORITHM 

— »T£NTER  atohs,  temperature,  etc,  ] 

■ A 

.VQMrWL  COMPOSITION  3Y  .MATERIAL  BALANCES 


^^2’  ^2* 

oxygen 

rich 

Ng,  Ar 

oxygen 

deficient 

COp,  CO,  H- 
Or  £ 
CO,,  H-,  H,0 

— ar  J 

c Z 2 

N,  IMATERIAL  BALANCE 


Hj  j C0M8INED  MOLE  NUM8ER  CONSTANT 

K2°  ' AND  MATERIAL  balance  in  hydrogen 


0.  H,  OH,  N,  so 
BY  MOLE  NUMBER  CONSTANTS 


oxygen 

rich 


CO  MOLE  NUMBER  CONSTANT  FOR  C02 
C02  MA-'IAL  BALANCE  IN  CARBON 
Oj  .MATERIAL  BALANCE  IN  OXYGEN 


oxygen 

deficient 


C02  MATERIAL  BALANCE  IN  OXYGEN 
CO  MATERIAL  BALANCE  IN  CARBON 
02  MOLE  NUMBER  CONSTANT  FOR  CO, 


RECYCLE 


( CHECK  CQR  CONVERGENCE^—-  . ■ 
j NO 

NEW  TRIAL  VALUE  FOR  02 
MEAN  OF  COMPUTED  AND  QLO  TRIAL  VALUES 


TO  PERCENT 
CALCULATIONS 
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Appendix  E 

INTERNAL  ENERGY  OF  FORMATION 


1.  The  heat  of  combustion  of  TN i is  listed  as  3620 
cal/g.  What  is  the  corresponding  value  tor  the  standard 
enthalpy  of  combustion  A£^  in  kilojoules  per  mole?  The  item 
ififa  includes  the  Greek  letter  A that  represents  a difference; 
subscript  a indicates  that  the  difference  is  for  a combus- 
tion reaction  between  the  fuel  and  oxygen.  The  letter  H in- 
dicates the  therraodyna.  c property,  the  enthalpy  (heat  con- 
tent). The  superscript  • indicates  a standard  value  with 
each  component  a pure  material  at  standard  temperature  and 
pressure — 25*C  and  1 standard  atmosphere. 

Answer^.  The  formula  for  TNT  is  CyHsNgQg,  and  the  cor- 
responding formula  mass  is  227.13  g/raol.  Then  on  a mole 
basis  the  above  heat  of  combustion  is  3620  x 227.13/1000  - 
822.2  kcal/mol.  In  tne  modern  thermodynamic  convention, 
where  negative  enthalpy  difference  indicates  a loss  of  en- 
ergy, Afl£  * -822.2  kcal/mol  * -3140.1  kJ/nol. 

2.  Presume  that  the  combustion  of  TNT  value  pertains 

to  a higher  heat  of  combustion  with  water  taken  as  a liquid, 
as  customary  for  unspecified  conditions,  and  find  a corres- 
ponding standard  enthalpy  of  formation  for  TNT. 

* 

Answer.  It  is  convenient  to  use  the  forn-'t  of  a chemi- 
cal equation  for  the  combustion  reaction.  Each  component  is 
then  assigned  its  standard  formation  value;  the  sum  for  the 
products,  Ho,  minus  that  for  the  reactants,  Hi,  becomes  the 
enthalpy  difference  for  the  reaction.  The  standard  forma- 
tion value  for  carbon  dioxide  is  -94.0513  kcal/mol,  that  for 
liquid  water  is  -68.3171,  and  the  values  for  oxygen  and  ni- 
trogen, of  course,  are  zero.  Writing  the  equation  for  the 
idealized  combustion, 

G7H5N3O6  + 5.25  02  - 7 CO2  + 2.50  H20  + 1 . 50  Nf3 

-94.0518  -68.3171 

A H°e  + 0 ■»  x7+  x 2 . 50  + 0 

J -658.363  -170.793 

' » >■ 

Hi  = A Hf 
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A - (H2  - » (-829 . ] 55  - iBj)  - -822.2  kcal/mol 

where 

IHf  - -829. 1S5  -(-822.2)  - -6.953  kcal/mol  - 

-29  100  J/ mol 

3.  Whai  issthe  corresponding;  standard  internal  energy 
of  formation,  bE*? 

4 

Anawar.  By  definition,  internal  energy  E is  the  en- 
thalpy  3 minus  the  pressure  .volume  product  PV . For  n moles 
of  gases,  this  product  is  nET . Thus  for  the  formation  re- 
action, the  internal  energy  of  formation  - LHf  ~ RTLn. 
For  TNT  formation, 

7 C(s)  + 2. 50  H2  + 1. 50  N2  + 3 02  - C7H5N306  (s) 

nj  - 2.50  + 1.50  k 3 * 7,  and  n2  ■ 0,  and  An  ■ -7 

Thus  hEf  • -29  100  + 8.314  34  x 23S.15  x 7 » - 11  750  J/raol . 

This  value,  along  with  others,  is  included  in  Table  E-l  of 
this  appendix.  Differences  indicated  there  are  not  unusual 
ones  and,  as  noted  above,  do  not  introduce  appreciable  un- 
certainty into  energy  computations. 
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TABLE  E-t.  Tabulated  A£,  Valuta  for  26SC. 


Acetone  

Acetone/nltrometltane,  90/M 

Acetylene  

Ammonia  

Baniana  

3 an  tana  

1 , 3-3u  tad  tana  


Carbon  dioxide 
Carbon  monoxxte 

Carbon  (a)  

Caiiuioaa  trinlirata 
Como.  3*3  

Cyanogen  

Cyclohexane  

Cydotol,  75/39  . . 
□iathyl  athar 
1.  2-Olnitrobeniene 
1,  3-Oinltrotoenxent 
1 , 4-Olnltrobentana 


Ethylene 

Ethylene  oxide  . 
Ethyl  nitrlta  . . . 
Glyceryl  trlnitrara 
(nitroglycerin) 

n-Haieana 

n-Haxana  '. . . . . 

Hexanitroben  tana 
Hexanittobentene 


Hydrogen  

Hydrogen  cyan  Ida 
Hydroxyl  


Mannitol  hamnltrata 

Mathana  

Mathanol  

a-Methyinechtheiana 
Monatomic  hydro  dan 
Monatomic  nitrogen 
Monatomic  oxygan 
Naohthaiana  


wm 

N 

0 

matt, 

{)/mo4 

6.00 

0 

1.00 

58.1 

4.93 

0.48 

1.48 

69.2 

2.00 

0 

0 

26,0 

3.00 

1.00 

0 

17.0 

6.00 

0 

0 

78.1 

6.00 

0 

0 

78.1 

6.00 

0 

0 

54.1 

1000 

0 

0 

58.1 

0 

0 

2.00 

44.0 

0 

0 

1.00 

28.0 

0 

0 

0 

12.0 

7.00 

3.00 

11.00 

297.2 

2.50 

2.19 

2.68 

100.0 

0 

2.00 

0 

52.0 

12.00 

0 

0 

84.2 

2.98 

2.39 

2.69 

100.1 

10.00 

0 

1.00 

74.1 

4,00 

2.00 

400 

198.1 

1.00 

2.00 

4.00 

168.1 

-.00 

2.00 

4.00 

188.1 

6.00 

8.00 

12.00 

464.3 

6.00 

2.00 

4.00 

182.7 

8.00 

0 

0 

30.1 

6.00 

0 

1.00 

48.1 

4.00 

0 

0 

28.1 

4.00 

0 

1 00 

44.1 

5.00 

1.00 

2.00 

75.1 

5.00 

3.00 

9 00 

227.1 

14.00 

0 

0 

86.2 

14.00 

0 

0 

862 

1 0 

6.00 

12.00 

348.1 

0 

8.00 

12.00 

34&1 

8.00 

8.00 

8.00 

296.2 

3.24 

2.02 

2.52 

100  0 

4.00 

2.00 

0 

32.1 

2.00 

0 

0 

2.0 

1.00 

1.00 

0 

27.0 

1.00 

0 

1.00 

17  0 

10.00 

0 

0 

118.2 

17.00 

0 

0 

125.3 

17.00 

0 

0 

113.3 

3.73 

1.69 

3.39 

100.0 

8.00 

6.00 

18.00 

462.2 

4.00 

0 

0 

16.1 

4.0 0 

0 

1.00 

32.1 

10.00 

0 

0 

142.2 

1.00 

0 

0 

1.0 

0 

1.00 

0 

14.0 

0 

0 

1.00 

16.0 

8.00 

0 

0 

128.2 

208  474 

• 148  751 
228  731 
-40  420 


79 


116  232 
393  510 
111  805 

0 

627  600 


2 


309  072 
141  362 


182  319 
181  478 


-9  5 


-402  771 

-400  000 
-104  140 
-608  966 


216  700 
358  780 
246  280 


■ vi-'1'- 
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TABLE  E-1.  (Contd.) 


Nun* 

Pormule 

Pdrmule 

meet, 

9/mol 

££*,  J/rhOl 

C 

H 

N 

0 

Nitric  odd*  . 

0 

0 

. 1,00 

1.00 

30.0 

90  370 

Nltroddruene  

9.00 

5.00 

1.00 

2,00 

123.1 

18  231 

Nhromeihene  

1.00 

3,03 

1,00 

2.00 

81.1 

-107  835 

NHromethiM/eutone,  78/26 

1.52 

3.7a 

0.74 

1.74 

100.0 

-1P8  929 

/rOctere  

8.00 

18.00 

0 

0 

114.3 

-1 38818 

Octol  

1.78 

2.58 

2.38 

2.99 

100.0 

20  210 

Oxv*en  . . • 

0 

0 

0 

2.00 

32.0 

0 

PBX-9007  

1.97 

3.22 

2.43 

2.44 

100.0 

3'.  859 

Ptnruvlvenie  cocl  du«t  

8.16 

5.26 

0.10 

0.55 

100.0 

-1C3  496 

Pen  ten*  

5.00 

12.00 

0 

0 

72.2 

-134  046 

Pemenltrobennene  

6.00 

1.00 

5.00 

10.00 

303.1 

80  000 

Penloli'e.  50/50  

6.16 

6.25 

3.41 

8.51 

334.4 

-97  903 

C/>rr'*rt4-Perhvdrofluo/ene  ... 

13.00 

22.00 

0 

0 

178.4 

-104  924 

petn  

5.00 

8 DO 

400 

12.00 

316.2 

-480  774 

Propen  • 

3.00 

8.00 

0 

0 

44.1 

-98  410 

ROX  

3.00 

6.00 

a oo 

6.00 

222.2 

111  471 

Rhode  Itlend  coel  duet  

6.80 

0.50 

0.01 

0.11 

100.0 

-90  749 

RJ4*J  

12.00 

20.00 

o 

o 

1643 

-113  585 

RJe&  

14  00 

18.40 

o 

0 

186.: 

41  235 

6.00 

10.00 

o 

5.00 

"82-2 

-948  371 

TACOT  

12.00 

4.00 

8.00 

8.00 

38M.2 

660  341 

TA  Z 

1 00 

0.00 

9.00 

o 

1472  | 

464  350 

Tetrenltrotoenxene  

6.00 

2.00 

4.00 

8.00 

258.1  | 

30  030 

Tetry!  

7.00 

8.00 

8.00 

8.00 

227.1 

38910 

wro-THOC  

10.00 

18.00 

0 

0 

138.3 

-104  924 

TNT  

7.00 

8.00 

3.00 

3.00 

227.2 

-11  750 

TNT  

7.00 

5.00 

3.00 

6.00 

227.2 

-13  877 

TNT  

7.00 

5.00 

3 00 

8.00 

227.2 

-12  170 

1,  2,  4-Trlttiwobeniene  

6.00 

3.00 

3.00 

5.00 

2131 

35  291 

1 . 3.  5-Trlnitrob«it»ne  

6.00 

3.00 

3.00 

8.00 

213.1 

-6  650 

a. oo 

o 

12.00 

6.00 

338.2 

1 180  000 

Weter  (II  

0 

2.00 

0 

1.00 

18.0 

-282  090 

Wetef  (vl  

0 

2.00 

0 

1.00 

18.0 

-240  590 
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Appendix  F 

THE  KOPP  RULE  FOR  ADDITIVE  ATOMIC  HEAT  CAPACITIES 
(CONDENSED  PHASES  NEAR  ROOM  TEMPERATURE) 


TABLE  F-l.  Additive  Atomic  Constants. 


16.7 

20.9 
22.6 

25.9 


! Atom 

J/ (mol • K) 

Atom 

C 

7.  5 

0 

I H 

9.6 

F 

1 5 

11.3 

P,  s 

; si 

15.9 

All  others 

The  heat  capacity  for  TNT,  formula  673^305,  is  esti- 
mated to  be  about 

(7  x 7.5)  + (5  x 9.6)  + (3  x 25.9)  + (6  x 16.7)  » 
278  J/ (mol • K) 
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Appendix  G 

INTERNAL  ENERGIES  FOR  COMBUSTION  PRODUCTS 

Standard  internal  energy  in  joules  per  mole  relative 
to  the  elements  at  25*C  as  a function  of  absolute  tempera- 
ture T in  kelvins.  Temperature  range  500-4000  K,  ideal  gases. 


£”*  - a + £(271000)  + o(T/1000)2  + i(  1000/2*) 


I 


I 

i 


TABLE  G-l . Coefficients  for  Calculating  Internal  Energies. 


Formula 

Coefficients 

a 

£ 

a 

d 

co2 

-423  070 

50  948 

589 

5280 

CO 

-126  080 

27  200 

327 

2440 

h2o 

| -264  870 

36  258 

1938 

5455 

h2 

-8  800 

20  898 

1327 

1112 

°2 

-11  085 

25  788 

994 

1176 

-14  275 

26  783 

370 

2508 

0 

245  030 

12  152 

76 

-252 

a 

213  050 

12  472 

0 

0 

0H 

28  533 

23  013  | 

929 

1725 

N 

468  340 

12  000 

0 

0 

NO 

76  195 

27  835 

260 

2247 

Ar 

-3  718 

12  472 

0 

0 

c(s) 

-7  500 

25  000 

0 

0 

I 

% 
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Appendix  H 

EMISSIVITY  COMPUTATIONS 


In  one  experiment  a 10.1%  mixture  of  methane  in  air  in 
a container  at  27*C  and  986  millibars  and  initially  quies- 
cent conditions  was  ignited  by  a central  electric  spark. 

The  blast  overpressure  generated  was  6.45  bars,  in  contrast 
with  the  adiabatic  value  computed  as  8.05  bars. 

1.  Compute  actual  products  temperature  and  internal 
energy  for  this  explosion  and  compare  with  those  for  adia- 
batic conditions. 

Anawer . By  iterative  calcuations  for  composition  as 
described  above,  and  for  series  of  trial  temperatures,  the 
temperature  corresponding  to  observed  overpressure  was  2161  K. 
The  corresponding  internal  energy,  relative  to  the  elements, 
at  25°C  is  -295  kJ/m3.  The  heat  flow  is  the  difference,  or 
801  kJ/m3  for  this  explosion. 

2.  In  centrally  initiated  internal  explosions  the 
flame  front  is  found  to  be  spherical  and  to  advance  at  a 
uniform  rate.11  For  the  above  explosion  the  rate  of  advance 
was  measured  as  0.333  m/s.  Compute  overall  flame  emissivity 
(the  ratio  of  actual  thermal  emission  to  that  of  a blackbody 
at  the  same  temperature). 

Answer.  The  Stef an-Boltzmanu  relation,  in  engineering 
format”!  is" 

Q - - *2>  (H-l ) 

where  Q is  the  rate  of  radiant  heat  transfer  and  a is  the 
Stefan-Boltzmann  constant,  5.67  x 10“8  W/m-K4.  Item  is 
an  area  factor  that  for  a completely  enclosed  flame  is  its 
emitting  surface  area.  For  a centrally  initiated  spherical 
flame  front  at  radial  distance  r from  point  of  initiation, 
the  emitting  area  is 

■ 4irr2 

Item  fg  of  Equation  H-l  is  a blackness  factor;  for  a flame 
completely  enclosed  by  cool  walls  this  is  gi  n as  the 
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